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Some effects of tfie anisotropy in a simple lattice gauge model at finite 
temperature 

L. A. Averchenkova, V. K. Petrov ^ 

'^Bogolyubov Institute for Theoretical Physics, 

National Academy of Sciences of Ukraine, Kiev 143, UKRAINE 

Monte Carlo simulations are carried out on the (3 + l)-dimensional Z(2) anisotropic lattice model, and a new 
method to simulate extremely anisotropic lattice systems with discrete symmetries is proposed. Dependence of 
the temporal and spatial average plaquette, Wilson loops on the anisotropy parameter is presented. 



1. Introduction 

In our previous paper 0], the phase struc- 
ture for Z(N) pure gauge theory with the lat- 
tice anisotropy ^ = aa/ar Q [aa-iar] is the spatial 
(temporal) spacing] has been studied analytically. 
It has been particularly indicated there that ar- 
bitrarily chosen trajectories in anisotropy param- 
eter and gauge coupling space can cross phase 
transition lines, and give rise to discontinuities in 
thermodynamic quantities along those trajecto- 
ries. 

The Z(N) gauge theory in four dimensions was 
first analyzed numerically for the Wilson action in 
a classical paper by Creutz, Jacobs and Rebbi 
Numerous Monte Carlo studies concern various 
aspects of Z(N) gauge theories, i.e. order of phase 
transitions both in pure gauge theories and 
in Z(N) gauge-Higgs models ||^-||], dependence 
of the phase picture on d |^ and on N 1^,^,^, 
elaboration of improved algorithms to simulate 
such systems with metastable states Our 
paper is aiined at clarifying the role of the lattice 
anisotropy ^ in the phase structure description of 
the finite temperature 4d Z(2) gauge theory by 
using Monte Carlo technique. Here, we check the 
analytic results of and picture the phase plane 
the area we cannot investigate analyt- 
ically. 



^ The bare anisotropy ^ = Kt/kct, KctCk-t) denoting the 
spatial (temporal) gauge couplings, is often considered as 
a proper lattice anisotropy parameter instead of g. 



2. Algorithm 

The partition function is a sum over all config- 
urations of the system 

Z = ^ exp(-Ki,^5); 



(1) 



s 



Re(CT; 



); 



= 0,---,iV-l. 



Although this is a finite sum, the number of con- 
figurations {<Jx;v\ is so large even for lattices 
which are only a few sites on a side that an eval- 
uation of the expectation value of a quantity O 



(2) 



cannot be performed directly. The Monte Carlo 
method replaces direct evaluation by generation 
of a sequence of configurations which simulates 
an ensemble of states in thermal equilibrium at 
inverse temperature Ki^^. 

We have used a heat bath algorithm which cre- 
ates a Markovian process as follows: a new value 
a'x-jj for the link variable is selected independently 
on the previous value of ax;u in a stochastic man- 
ner with the probability 

y^{qx-u -> q'x-u) ^ exp{-K^f,So{q'x-^)), (3) 



(4) 



2 



(5) 



exp {27ri(qx+,y.^ - qx+i_i-u - qx-,i)/N) . 



The transition matrix W{qx;v I'x u) obeys an 
equation of detailed balance 



W{qx;u -> q'x;u) ^ {-«:„^(5o(9L;J-5ofe;.))} 



(6) 



which is a sufficient condition for final distribu- 
tion to be the Boltzmann one. A complete cycle 
through all the link variables of the lattice when 
the sites are chosen in a some way (in our case, in 
a random way) is called as a Monte Carlo sweep. 
The detailed balance (||) leads to that the prob- 
ability to find a configuration {cr} after n ^ oo 
sweeps 



For the Z(2) gauge group 



'PWx-u = '^) = c-exp - ^ Ki,^a^.^^ 

"^(f^i;!/ = ~1) = C-eXp ^ Kj,^(Tj;;^p 

Since the Z(2) gauge group is a discrete group, 
-6 < 'S'x-iyt^ < 6, the sums Y^^ax-nm and 

^x-riT take the finite number of values, thereby 
can play the role of an index in some matrix 



^„,s, = exp 



C^er ^ ^ ^x;nm ^ ^ <^x\nT ] (7) 

cr T j 



Via' = 



(8) 
(9) 



Neither exponents nor products (expensive com- 
puting operations) are calculated in the course of 
sweep, they are computed beforehand. Actually 
we calculate only the sums 'Ylia ^'^'^ ^"^^ 
out the corresponding preliminary calculated ex- 
ponent in the table (|^) as the probability to find 



a configuration {cr}. This algorithm is applied 
to every link of the lattice, giving us a next in 
turn configuration from the Boltzmann distribu- 
tion set. This procedure is repeated many times, 
and an estimation of (O) is obtained. 

There are, however, many problems - both 
standard and specific for the discrete groups - 
in applying the above algorithm. Due to a dis- 
cretity of the Z(2) gauge group, metastable states 
appear, and the system will then remain there for 
many sweeps if the tunneling probability between 
states is small, especially at large Ko-^r- This prob- 
lem can be hardly solved by increasing the num- 
ber of sweeps, and in some area of parameters a 
new algorithm is required. We have also written 
the procedure where new and previous states dif- 
fer in the value of all link variables which adjoin 
the selected link ax\v The correlation between 
successive configurations then is supposed to be 
weaker. Here, we present the results obtained 
by using a standard heat bath algorithm. Any 
procedure gives (0) with statistical errors. The 
measurements cannot be considered as indepen- 
dent, of course, because they are obtained on the 
basis of a sequence of highly correlated config- 
urations. By taking into account a strong cor- 
relation between the measurements (which is an 
inevitable consequence of a "microscopic" nature 
of the algorithm - one link changes at one time), 
the statistical error is computed as 



= Vx2(l+r)/(l~r), 

= ^(O,-0)V(ri-l)n, 

i 



(10) 



where n is the number of samples in the aver- 
age. For discrete groups, especially at large Kjy^, 
the problem of correlations becomes to be very 
sharp, because we can stick at one configuration, 
the correlation becomes then practically infinite 
and the measured average is questionable. 

3. Results 

We have performed the simulations on the lat- 
tice of N1 X Nr (12^ X 4) size, the total num- 
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ber of sweeps is equal to 500000 x 2, we start 
from a completely disordered configuration, two 
thousands sweeps were used for thermalization. 
Various quantities such as the space-like (time- 
like) average plaquette {Pc,){{Pt)), Wilson loops 
W{I,J)a,T, the Creutz ratio x{I,J)[I ^ J ^ 
2, 3, 4] have been estimated as functions of /3 = 
y/KaKr = 4:/ g'^ and 

The temporal W{2, 2)t (crosses in Figs) and 
spatial W{2, 2)cr (circles in Figs) Wilson loops at 
different ^ are plotted here. When increasing /3, 
the Wilson loop becomes nonzero first in the tem- 
poral direction and next in the spatial directions 
- which is not surprising because ^ enhances the 
time-like plaquettes and suppresses the space- like 
plaquettes. As is well-known, at high temper- 
ature (this corresponds to the high anisotropy 
f ^ 1) the time-like Wilson lines acquire nonzero 
expectation values, while the space-like Wilson 
lines do not. 

Data for the Wilson loops {Wa,T (I x J)) ~ 

exp{— acr,T (/^'^) ^ J)} '^^^ be roughly fitted 
by 

a,,r (/?, - 2 (1 + £.,r) • - /3) e - /?)} 

with < e^^r < 1; P^^r = 6^"".-; ^ 0.5, Vr ^ 
—0.5 -which does not agree with the condition 
obtained for Wcr,r in [|| for SU(N) gauge group. 
In particular, W^.t cannot be fitted by the uni- 
versal function / (^/3, ^ = = C ■ 
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Figure 1: C = 0.05 
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Figure 2: ^ = 1 
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Figure 3: ^ = 10 
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